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Abstract 

A calculational approach in fluid turbulence is presented. Use is 
made of the attracting nature of the fluid-dynamic dynamical system. 
An approach is offered that effectively propagates the statistics in 
time. Loss of sensitivity to an initial probability density functional 
and generation of stationary statistical effects is speculated. 

1 Introduction 

Turbulence of a system is thought to be a condition of the system for which 
the space and time evolution of the dynamical variables are chaotic, but for 
which a statistical distribution functional is applicable. The time evolution 
of an initially given probability distribution functional for the dynamical 
variables is given by the deterministic evolution equations that govern the 
system, called deterministic chaos. 

A fluid environment is a dynamical system, having energy input mech- 
anisms at long distance scales, and energy dissipation mechanisms due to 
viscosity occurring at small distance scales. Hence, the fluid has an attrac- 
tor, called the Strange Attractor in the turbulent regime, for its dynamical 
variables. These considerations apply to the standard system of interacting, 
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or self-interacting, classical fields which are governed by evolution equations 

that are at most first order in the time derivatives, and that possess dissi- 
pation mechanisms, together with stationary boundary conditions. It is pro- 
posed here that the fluid probability density functional also has an attractor 
for its time evolution, and an approach to generating this time evolution is 
presented. The same mechanism that causes the dynamical variables to have 
an attractor in phase space, that is, the tendency for the equilibration of 
energy input rates and energy output rates to set in, also causes an arbitrary 
initial statistics to evolve toward an attracting statistics, which is stationary 
in time. 

It is the stationary statistics that allow the Kolmogorov scahng ideas to 
have applicability. The evolution of the fluid's statistics can be set up as part 
of a space-time path integral. Ensemble averages of any dynamical variable 
can be formulated in terms of this path integral. Fluid space-time sampling 
techniques naturally suggest a useful way, using a relatively arbitrary initial 
statistics functional, to calculate averages. 



2 Description of the Mathematical Approach 
to Turbulence 

Let us set up the evolution equations for incompressible fluid dynamics. The 
extension of the following to compressible fluids will pose no unusual diffi- 
culty. 

We have, 

dv -> 2 ^ 
P-^ = J + riV V, 

where / is an external force density, such as due to a scalar pressure field, 
and ?7 is the coefficient of viscosity, v is the fluid velocity field. We also have 
from the conservation of mass. 

Here, p is the mass density of the fluid. If this mass density is a constant, 
then the velocity field is divergenceless. Also, 

dv dv ^ „^ 
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So we have the fluid dynamic system, 



dv 
di 
V ■ V 







VP 



p 



(1) 
(2) 



where P is the pressure, and v = ^. 

We drop the external force density in what follows. This is not an essential 
step. What are needed are a set of interacting fields, together with stationary 
boundary conditions to allow a deterministic time evolution of the set. We 
also associate with the spatial velocity field a probability density functional, 
p[v,t\. The fluid statistics time-evolves according to deterministic chaos [1] 



That is, the number, p, associated with the spatial velocity field vq at time 
to, will be associated with f/ at time t/, where f/ is the velocity field vq 
deterministically evolves into from time to- Given a functional of the spatial 
velocity fields, A[v\^ its ensemble average, at time tf is. 



We want to propagate the fluid's statistics according to deterministic 
chaos, even though the detailed fluid orbits are chaotic. Let, 



0, 





where the kernel is the delta functional. 



K[vf,tf]Vo,to] = 5[vf - f[tf;vo,to]]. 



< A[v] > 




< A[v] > 




where. 



K[vf,vo] 



^[vf - flM] 

1 5b/-/iN]^K-/iN]rfN 



3 



= y S[Vf - f2[vi]]S[vi - f2[v2]]S[v2 - f2[vo]]d[vi]d[v2] 

= J 5[vf - /sbill^K - f3[v2\]S[v2 - /aball^bs - /sNjc^Klc^Nc^ba] 
' ' ' ) 

where the velocity fields, vi, V2, vs, ... occur in chronological order, vi being 

closest to the time tf. Eventually, we have an Jm, where M is large, such 
that vm = /m[i^o] is infinitesimally difi^erent from vq. 
Hence, 

<A[v] >= I d[v]A[vf]6[v - fM[v]]p[vo,to], 
where the functional integration measure is, 

d[v] = dfti/ldf^;!] • • • d[^;jM-]rf[t'o], 

and 

< A[v] >= J A[fM[- ■ ■ [/m[/mN]] ■ ■ MvoMdH- 

The exact rule, /m^], requires a solution of the fluid dynamic equations, 
incorporating the boundary conditions. The exact rule, JmIv], is difficult 
to find. Let us use an approximate rule, Fm[v]. Then, we may say, with 
motivation to follow, 

< A[v] > = / d[v]A[vf]5[v - Fm[v]]X[V ■ v]X[v - vb]p[vo, to]. (3) 

The functional integration is over all space-time velocity fields within the 
spatial system, between times to and tf. 6[v — Fm[v]] is a space-time delta 
functional. FMif] generates v from a an infinitesimal instant earlier. The 
A functionals are evaluated at a particular instant. The needed properties of 
the A functionals are A[0] — 1, and X[g ^ 0] — 0. X[g] could be, 

lim e-^'/'i. 

ei-*-0+ 

We have, then, 

< A[v] > = J d[vo]A[FM[FM[- ■ ■ [Fm[vo]] ■ ■ "IIIpK, ^o] 

X[Fm[- ■ ■ [Fm[vo]] ■■■]-vb]--- X[vo - vb] 

A[V • Fm[- ■ ■ [Fm[vo]] •••]]••• A[V • Fm[vo]]X[V ■ vo]. (4) 
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The right-hand side of (jlj) equals < A[v] >, because p[fo,to] "will be non-zero 
only for spatial fields Vq that make all the A's equal one. This means that we 
only utilize divergenceless spatial velocity fields satisfying, also, the spatial 
boundary conditions. 

The spatial velocity fields have an attractor, determined by the stationary 
boundary conditions on the fluid [3]. When the boundary conditions allow 
steady laminar flow to become established, the attractor consists of a single 
spatial velocity field. When the Reynolds number becomes large enough, 
bifurcations set in, or the onset of instability occurs, and the attractor begins 
to consist of more than one spatial velocity field. In the turbulent regime, 
the attractor consists of many spatial velocity fields, and the fluid accesses 
these fields according to a probability distribution |1] [5] [H]. 

Given a functional of the spatial velocity fields, A[v], and the fluid dy- 
namic system of equations, ([1]), (|2]), we will say that its ensemble average 
when the system has reached its attractor is, 

< A[v] > = lim / d\v]A[vf]S[v - FmM]A[V ■ v]X[v - vb] 

tf-to^CX> J 

■pbo,^o]- (5) 

The delta functional condition, — FmH] implements equation ([T]). Fm[v] 
is an approximate rule for carrying v from, an earlier instant to a later instant. 
In the first approximation, 

Fm[v] = V - At{v ■ Vv - uV'^v), 

where v is evaluated at an instant At earlier. -FM[f] can be expressed in 
all higher powers of At, with only higher order spatial derivatives being 
required. This is due to the fact that the hydrodynamic evolution equations 
are, at most, first-order in the time derivatives. A[V ■ f] implements a zero 
divergence condition on the spatial velocity fields, and X[v — vb] requires the 
spatial velocity fields to have values on the spatial boundaries. 

Let us consider the path integral ([5]) to be on a space-time lattice. We 
could use, 

S(x) = lim — 

We have for the average of A[v] in the steady-state (attractor), letting ei(Ax)^ = 
e, where ei is in the A-functional for the zero divergence condition, and Ax 
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being the Cartesian spatial distance step size, 



< A[v] > = lim lim / d[v]e-^^''^^'A[vf]p[vo]. (6) 

e-*-0+ tf-to^oo J 

Also, 

p[vo\ = pK,^o], 
and the functional integration measure, d[v], is, 

H[v] is a functional of the lattice space-time velocity field. Also, neglecting 
boundary effects, 

H[v] = Yj(('"i ^ ^"i-^ ~ 9[vi-i]^tf + (v^^ijki - v.^^i_ijki ^ Y) 



where 



or, 



g\v^ = -vi ■ Vvi + uV'^vi, 



r 1 {Vijkl — Vi-ljkl) , , (Vijkl — 2Vi-ijkl + Vi-2,jkl) , 
g\Vi\ = -Va:inkl 1 ':!—L^ yv^^ Tir^ ':!—L^ . 

N is the number of spatial lattice points in the system, and T is the number 
of time slices. We have as a sum over all space-time lattice points in the 
system, and ^ a.s a. sum over all space-time lattice points on the spatial 
boundaries. Also, ijk are spatial lattice indices, and I is the index in the 
time direction. 

This discretization technique is expected to get better as one increases 
the lattice fineness and makes use of higher powers of the time step and 
higher order finite difference approximations of the spatial derivatives 
A good approximation to the attracting statistics as a starting point will 
shorten the evolution time required for averages occurring in the steady state. 
Gaussian statistics, however, could be a good generic starting point. The 
path integral ([6]) can be evaluated with Monte Carlo techniques utilizing 
importance sampling. A calculation of the stationary spatial velocity field 
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that would exist, for the given boundary conditions, if that field were stable, 
is expected to be a good starting point from which to begin a sampling of 
the velocity field space-time configuration space. The average values, for the 
starting Gaussian statistics, can be taken as the values of the velocities of 
the calculated stationary spatial velocity field. 
We have, 

<A[v] :,^j±^hApM^ 

where, 

and the space-time configurations are distributed according to the weighting 
g--f^H/<:_ is the number of space-time configurations in the sample. 
is the observable A[v] evaluated at the final time slice of the i*^^ space-time 
configuration. The value p[t'o,i] is attached to the initial time slice of the i*'* 
configuration. We need. 



So, we must do the integral, J, and constrain J2 p['Vo,i] to equal n/ J, possibly 
by varying the variances in the Gaussians of p[v, to]. 



3 Summary 

We have said that the time evolution of the statistics also has an attractor 
for the following reasons; (1) It is a way to get a solution to the problem of 
arriving at the steady state turbulence statistics. One knows that the steady 
state statistics is stationary with respect to its time-evolver. Postulating that 
this statistics is the attractor of its time-evolver means one does not have to 
have the statistics to get the statistics, thereby offering an approach to the 
closure problem for the determination of correlations. 

(2) The statistical physics approach has been successful in equilibrium 
thermodynamics, where the derivation of the microcannonical ensemble can 
be taken as the indication that equilibrium is the attractor for the dynamical 
system when there is a Galilean frame for which the boundary conditions 
allow no energy input into the system. In the attractor, the mean energy 
input is equilibrated with the mean energy output, because in the attractor 
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dissipative losses have effectively shut off, and the system becomes effectively 
conservative. That is, the system becomes describable by a time-independent 
Hamiltonian. The stationarity of the statistics requires the vanishing of the 
Poisson bracket of the statistics with this Hamiltonian, resulting in the statis- 
tics of equal a priori probability. 

(3) The dynamical system, of which a fluid is an example, has an attractor 
[3]. The dynamics of the statistical approach should mirror the dynamics of 
the actual dynamical system. 

(4) The statistics of the dynamical system, prior to reaching the attractor, 
has no reason to be unique. The statistics of the attractor is expected to 
be unique, in which the geometry of the system, the stationary boundary 
conditions, and viscosities, all of which determine the Reynolds number, play 
a crucial role in determining the attractor. 

(5) When the stationary statistics of the fluid is reached, the equilibration 
of energy input and energy output rates has set in [8]. 

4 Conclusions 

In the discretized version of the path integral that attempts to arrive at the 
stationary statistical effects in the generation of the ensemble average of an 
observable, using an approximate rule for the dynamics, one should arrive at 
a greater insensitivity to a starting statistics, and a generation of stationary 
statistical effects, to within any desired level of approximation. One is trying 
to use the turbulent transience to get at steady state turbulence effects. 
These stationary statistical effects become the backdrop for Kolmogorov's 
ideas of self-similarity and the resulting scaling exponents. 

The potential usefulness of lower dimensional models should not be un- 
derestimated |9]. However, it is hoped for, here, that a transition to a detailed 
use of the Navier-Stokes equations can be brought about. 
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